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Abstract. The detection of TeV 7-rays from the blazar H 1426-1-428 at an integral flux level of (4 ± 2stat ± Isyst) x 
10~^'^ergcm~^s~^ above 1 TeV with the HEGRA imaging atmospheric Cherenkov telescope system is reported. 
H 1426-1-428 is located at a redshift of 2=0.129, which makes it the most distant source detected in TeV 7-rays 
so far. The TeV radiation is expected to be strongly absorbed by the diffuse extragalactic background radiation 
(DEBRA). The observed energy spectrum of TeV photons is in good agreement with an intrinsic power law 
spectrum of the source oc E~^'^ corrected for DEBRA absorption. Statistical errors as well as uncertainties about 
the intrinsic source spectrum, however, do not permit strong statements about the density of the DEBRA infrared 
photon field. 
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1. Introduction 

Many nonthermal extragalactic objects representing dif- 
ferent classes of AGNs are considered as potential sources 
of TeV photons. First of all this concerns the BLLac 
population of blazars in general; two nearby representa- 
tives of this class, Mkn421 and MknSOl with redshifts 
of 2;=0.030 and z=0.034, respectively, are firmly estab- 
lished as TeV 7-ray emitters. Of special interest are the 
so-called "extreme synchrotron blazars", BLLac objects 
with flat spectra of synchrotron emission and high X-ray 



to radio flux ratios (e.g. Costamante et al. 2001). The ob- 
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servations of two such objects, Mkn501 (Plan et al. 1998) 
and lES 2344-^514 ( piommi et al. 2000] ), by BeppoSAX 
showed that in a flaring state the synchrotron peak of 
"extreme blazars" can reach 100 keV. Remarkably, both 
these objects are also repo rted as TeV blazars (see e.g. 
[Catanese fc Weekes 1999 ). Most probably, this is not 
a mere coincidence; since the high synchrotron peak is 
an indicator of acceleration of electrons to ultrarelativis- 
tic energies, the "extreme blazars" are obvious candi- 
dates for TeV emission. Consequently, some of these ob- 
jects were intensively monitored with imaging atmospheric 
Cherenkov telescopes (lACT) over the last several years. 
Possible detections of TeV signals have been claimed by 



different groups for lES 2344+514 (Catanese et al. 1998) 
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lES 1959+650 ([Nishiyama et al. 1999| ), and PKS 2155-304 
(Chadwick ct al. 1999). Most recently, statistically signifi- 
cant TeV signals have been reported also for H 1426+428, 
by the VERI TAS (|Horan fc the VERITAS coll. 2001D 
and HEGRA ( potting et al. 200l| ; [Aharonian 200l| ) col- 
laborations. The relatively large redshift (z=0.129) of 
H 1426+428, as well as the fact that it is the third 100 



: angular cut 



keV blazar ever found (Costamante et al. 2001), make 
this object an extremely important target for future mul- 
tiwavelength studies of the nonthermal processes of parti- 
cle acceleration and radiation in AGN jets, as well as for 
probing the DEBRA in the near infrared region. 

In this letter we report the HEGRA observations of 
H 1426+428 performed in 1999 and 2000, and briefly dis- 
cuss some apparent astrophysical implications suggested 
by the obtained results. 

2. Observations and results of analyses 

For the following analysis, 14.5 hrs of data from 1999 and 
29.9 hrs from 2000 were available, taken with the HEGRA 
lACT system; see Table |l| for the dates of observations. 
The median zenith angle was 17°, resulting in a 7-shower 
peak detection rate at an energy of 700 GeV for Crab-like 
spectra (Konopelko et al. 1999). 

Shower reconstruction and cuts were described e.g. in 



Aharonian et al. (1999a). The methods applied in this 
analysis differ slightly, depending on the tasks which were 
performed. For the search for new sources we applied tight 
cuts and stereo algorithm #3 ( [Hofmann et al. 1999 ) to 
achieve high sensitivity. For the spectral analysis of the de- 
tected source, somewhat relaxed cuts were used which pro- 
vide a larger 7-acceptance; these cuts are less susceptible 
to systematic errors, and the spectral reconstruction was 
tested in great detail using Mo nte Carlo studies and 7- ray 
data from known TeV sources (lAharonian et al. 1999a| ,|b|). 



1999: Feb. 23., Mar. 15., 16., 18.-22., Apr. 10.-12., 20. 
2000: Mar. 5.-7., 14., 28.-30., Apr. 3.-5., 26.-28., 30., 
May 1., 6.-8., 29., 30., Jun. 1.-5. 



Table 1. Dates of observations. Each date refers to the 
respective following night, with typically one hour of ob- 
servations. 
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Fig. 1. Dots: Number of events vs. the squared angular 
distance to the position of H 1426+428. Shaded histogram: 
Background estimate; up to 9^ — 0.0225 deg'^, data from 7 
control regions are used. Therefore, the statistical error of 
the background estimate is much smaller than the error of 
the source distribution. The dashed line shows Crab excess 
events, measured at similar zenith angles, scaled down to 
6%, and superimposed on a flat background. The vertical 
dotted line indicates the position of the optimum angular 
cut, based on Crab data. 
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Fig. 2. The sky area around H 1426+428 in celestial co- 
ordinates (epoch 2000). The circles show the 1 and 2 a 
error of the reconstructed TeV position (HEGRA data). 
The systematic pointing uncertainty of the HEGRA sys- 
tem is 25 arcsec. The positions of X-ray selected QSOs (t^: 
H 1426+428 (position from Hubble space telescope obser- 
vations), A: CRSS J1428. 3+4231, T: CRSS J1429. 7+4240) 
and galaxies (•: the galaxy group CRSS J1429. 1+4241) 
were obtained from NED. 





signal 


spectrum 


stereo algorithm 


#3 


#1 


shape cut: mean scaled width < 


1.1 


1.2 


angular cut: < 


0.011 deg^ 


0.05 deg^ 


iVoff 


1779 


6258 


a = Aon/Aoft 


0.143 


0.262 


Non 


360 


1839 


A^ 7 — candidates 


105.9 


199.2 


significance a 


5.8 


4.3 



Table 2. Cuts, event numbers, and significances for the 
signal search and the spectral analysis. A: on resp. off 
source area. 



The background was estimated using a set of control re- 
gions in the field of view (FOV) ; the setup provides seven 
times more background data th an on source data unde r 
identical acceptance conditions ( Aharonian et al. 2001 ). 
The optimum angular cut (for the signal search) was de- 
rived on the basis of the angular resolution - using 7-ray 
events from the Crab nebula - and the background level. 
Table |^ summarizes cut parameters, resulting event num- 
bers, and significances for the signal search and spectral 
analysis. 

Figure |l| shows the event distribution, obtained from 
the signal search, as a function of the squared angular 
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Fig. 3. Differential energy spectrum of H 1426+428 (dots); 
the numbers indicate the (background subtracted) event 
counts and their corresponding statistical Icr-errors. A 
power law fit to the data is shown by the dashed-dotted 
line. The solid and the dashed lines represent absorbed 
model spectra, obtained by an initial power law spectrum 
oc E~^'^ (normalized to the data) combined with different 
absorption models (see Fig. 0) . 



distance to the source position. The excess significance, 
calculated according to Li & Ma (1983), amounts to 5.8 cr 



(2.4 cr in 1999, 5.3cr in 2000). While there is an indication 
that the source was in a higher state in May 2000, the 
total data set is also statistically consistent with constant 
signal accumulation. 

The event distribution in the FOV was used to re- 
construct the source position as seen in TeV 7-rays; 
for details of the procedure see Piihlhofcr et al. (1997 ) 
and Aharonian et al. (2001 ) . Figure ^ shows the sky 
area around H 1426+428. The positions of X-ray selected 
QSOs and galaxies were obtained from the NASA/IPAC 
Extragalactic Database (NED). The median angular res- 
olution of the HEGRA system is approx. 0?1; the given 
statistics allowed the reconstruction of the center of TeV 
emission with a statistical error of 50"(lcr). A confusion 
of the TeV source with other known X-ray sources is ex- 
cluded. We note that the very prominent X-ray source 
GB 1428+4217 is located at an angular distance of 41.'2 
to H 1426+428, far outside the area shown in Fig. ^ 

The energy reconstruction, which has a single event 
resolution of AE/E = 20%, and the spectral evaluation 
were described in detail in Aharonian et al. (19994 [fc|) ■ The 
energy spectrum was derived from the raw, background 
subtracted photon count spectrum using a so called effec- 
tive area, which is a response function depending upon the 
reconstructed energy and zenith angle; the effective area 
was adjusted regularly according to the varying detector 



conditions and different system setups (3/4/5 telescopes 
included in the system). 

Figure || shows the reconstructed differen- 
tial energy spectrum of H 1426+428. To first or- 
der, the spectrum can be described by a pure 
power law d^/dE = (2.0 ± l.Bstat ± O.lgyst) x 
10-i2(£;/TeV)(-2-6±0-6.t.t±o.w.t) phcm-2 s^^ TeV"^ 

(the systematic errors do not include the 15% error on 
the energy scale). The spectral slope is similar to the one 
measured for the Crab nebula (Aharonian et al. 2000). 
The measured integral flux (derived using loose cuts) 
above ITeV is (1.7±0.5stat iO.l^yst) x lO^i^phcm-^ g-i 
which corresponds to 10% of the Crab nebula flux. A 
comparison of the total photon count rates of H 1426+428 
and Crab nebula after tight cuts yields 6%; within 
statistical errors, both results are compatible. 

3. Discussion 

For any reasonable DEBRA model (see Fig. ^a), the TeV 
radiation from H 1426+428 is expected to arrive with a 
drastically modified spectrum. Reasons are the strong en- 
ergy dependence of the mean free path of 7-rays in the 
intergalactic medium, A{E), and the fact that the opti- 
cal depth r = cz Hq^ A~^{E) significantly exceeds unity 
at all 7-ray energies above 300 GeV for H 1426+428 with 
a redshift of z=0.129. The conclusion about strong ab- 
sorption at all energies covered by HEGRA does not de- 
pend on the large systematic uncertainties in the reported 
DEBRA fluxes. This is demonstrated in Fig. ^b where 
we show the intergalactic absorption factors calculated for 
two reference DEBRA models. Model 1 is close to the re- 



cent calculations of [Primack et al. (2001| ) based on the 
so-called Kennicut initial mass function (IMF), and fits 



quite well the DEB RA fluxes reported recently by [Wright 
|fc Johnson (2001] ) at 1.25, 2.2 and 3.5 /xm wavelengths. 
Model 2 is designe d to match the NIR fluxes reported by 
Matsumoto (2000 ) and Cambresy et al. (2001 ) who claim 
very high fluxes at wavelengths shorter than 2 /im. Despite 
significant quantitative differences, the two DEBRA mod- 
els result in similar absorption features in the TeV spec- 
trum: both models predict a very strong steepening of the 
spectrum at energies below 1-2 TeV, but an almost en- 
ergy independent absorption at energies around several 
TeV (see Fig. ^b). This effect is common for all realis- 
tic model curves describing the DEBRA in the near IR 
(e.g. Primack et al. 2001), since they reflect the character- 
istic shape of the starlight spectrum which in the wave- 
length band between 1 and several microns behaves as 
vl{v) oc A'' with /3 ~ — 1. The corresponding number den- 
sity of background photons n(e) oc results in a nearly 
constant optical depth for 7-rays at energies around sev- 



eral TeV ( [Aharonian et al. 1999bD . 

The X-ray spectrum of H 1426+428 which extends up 
to 100 keV favours an intrinsic pure power law spectrum 
in the relevant energy range, with an index of 1.9 (a de- 
tailed discussion is beyond the scope of this paper and will 
be given elsewhere). The observed spectrum - after ab- 
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Fig. 4. (a) Diffuse extragalactic background radiation. 
Tlie reported fluxes at far infrared wavelengths above 
60 ^m, as well as at optical wavelengths below 1 /im are 
taken from the recent review by Hauser fc Dwek (2001| ). 



The lower flux limit at 15 fim is derived from the ISOCAM 
source counts ( pranceschini et al. 2001 ). The fluxes re 
ported at NIR are shown by filled squares (Wright & 



Johnscn 2001), open circles (Matsumoto 2000), and open 
diamonds (Cambresy et al. 2001). Two model curves are 



shown by dashed and solid lines, respectively, (b) The 
spectrum modification factors exp[— r(i?)] calculated for 
the two DEBRA models for Mkn 501 and H 1426+428. 

sorption - is therefore not expected to obey a power law. 
The results presented in Fig. || show that the HEGRA 
spectral points are indeed better described by an intrinsic 
power law spectrum modified by intergalactic absorption. 
However, a power law fit with an adjustable spectral in- 
dex results in a reduced x^/dof = 4.4/3 = 1.47, while 
absorbed spectra - with a fixed source spectrum oc E~^-^ 
- yield xVdof = 3.7/4 = 0.9 and xVdof = 1.0/4 = 0.25, 
for the DEBRA models 1 and 2 (see Fig. respectively. 
Hence, more accurate spectral measurements, as well as 
detailed theoretical studies (concerning both the forma- 
tion of the 7-ray spectra in the source and their subsequent 
deformation during the passage through the extragalactic 
photon fields), are required. 

Nevertheless, the sheer fact of detection of TeV 7-rays 
from H 1426+428 should lead to a revision of the cur- 
rent conceptual view of TeV blazars, according to which 
the synchrotron (X-ray) peak in the SED dominates over 



the inverse Compton (TcV) peak (see e.g. Fossati et al 
1998). Indeed, although the detected energy flux of 7-rays 
is only ~ 4x 10^^^ crgcm^^s^^, corrected for intergalactic 
absorption this flux may well exceed 10~^''ergcm~^s~^. 
For comparison, the X-ray flux measured by BeppoSAX 
(postamante et al. 200 1|) and ASCA (T. Takahashi, pri- 



vate communication) is well below 10"^" ergcm~^s~^. 
Since the corrected TeV luminosity seems to exceed the 
level anticipated from the current models of TeV blazars 
by far, this result may have a crucial impact on the further 
development of models for TeV blazars. 

H 1426+428 is the most distant source detected in TeV 
7-rays so far. The HEGRA collaboration intends to con- 
tinue observing this source extensively in 2002 in the hope 
to decisively improve the statistics of the measurement. 
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